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ABSTRACT 

We demonstrate observational evidenee for the oeeurrenee of eonveetively driven 
internal gravity waves (IGW) in young massive 0-type stars observed with high- 
preeision CoRoT spaee photometry. This evidenee results from a eomparison between 
veloeity speetra based on 2D hydrodynamieal simulations of IGW in a differentially- 
rotating massive star and the observed spectra.We also show that the velocity spectra 
caused by IGW may lead to detectable line-profile variability and explain the occur¬ 
rence of macroturbulence in the observed line profiles of OB stars. Our findings pro¬ 
vide predictions that can readily be tested by including a sample of bright slowly and 
rapidly rotating OB-type stars in the scientific programme of the K2 mission accom¬ 
panied by high-precision spectroscopy and their confrontation with multi-dimensional 
hydrodynamic simulations of IGW for various masses and ages. 

Subject headings: Asteroseismology - Line: profiles - Stars: massive-Stars: rotation 
- Stars: oscillations (including pulsations) - Stars: waves -Techniques: photometry 


1. Introduction 


The existence of high-order gravity-mode (g-mode) oscillations in Slowly Pulsating B (SPB) 
stars, which are core- hydrogen burning stars with mass between roughly 3 and VM© (SPBs here¬ 
after, [^elkens|[l99lj) was established more than two decades ago, prior to understanding their ex¬ 
citation mechanism. We now know that these oscillations are driven by the he at-mechanism associ- _ 

ated with the opacity bump due to iron-group elements in the stellar envelope (IDziembowski, Moskalik, & Pamyati 
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1993l:lGautschy & Saidll993h . However, the detection o£g-mode period spacings associated with 


those standing waves as predicted from theory ( Tassoullll980 ~) and required for asteroseismology of 
such stars, r emained impossibb from the ground, eve n after carefully planned long-term dedicated 


campaigns (lAerts et al. 


19991: iDe Cat & AertsI1200211 . This is partly due to the low amplitudes of 


g-modes and is also made difficult because the periods of gravity modes in massive stars are of or¬ 
der days and these timescales have strong contamination by daily aliases in the amplitude spectra 
of ground-based photometry and high-resolution spectroscopy. Gravity-mode asteroseismology of 
SPBs only saw its birth due to the uninterrupted high-precision space photometry assembled with 
the CoRoT and Kepler missions. These data led to the detection of period spacings caused by 
dipole mod es of consecutiye rad i al order and offered the mode id entification necessary for seismic 
modelling ( Degroote et al. 2010 : Papics et al. 2012 . 2014 . 2015 ). 


While heat-driven g-mode oscillations in SPBs are strictly periodic and have a well-known 
and quantifiable effect on observed time -series photometry and spectral line-profile variations (e.g.. 


De Cat & Aerts 

20 

02; 

Aerts et al. 

2014 

), the effect of waves excited stochastically by core convec 

tion ( 

Samadi et al. 

bolO; 

Belkacem. Dupret. & Noels 

2010: 

Shiode et al. 

2013 

) on such observa 


tions is relatively unknown. This lack of observational diagnostics connected with IGW is par¬ 
ticularly relevant in the context of the variability of 0-type stars, in which heat-driven modes are 
not excited, while their large convective cores likely drive IGW efficiently. This work is therefore 
focused on the search for IGW driven by core convection in photometric and spectroscopic obser¬ 
vations of a few carefully selected observed 0-type stars. Such signatures are important because 
the existence of these waves could point to enhanced angular momentum transport and chemical 
mixing and hence, guide inclusion of these processes into future theoretical models. 


2. Modeling convectively driven waves 


There have been relatively few theoretical predictions for th e spectra and amplitud es of IGW 
excited by core convection in massive stars. On the one hand. Browning et al.l (120041) provided 
the first 3D simulations of the core convection for a 2 Mq star but these only covered the inner 
30% in radius and hence, omitted much of the wave propagation region. More recent theoretical 
models do consider the entire radius of the star but are limited to one-dimension (ID), neglect 
rotation and must make assumptions a bout the nature of convection and convective-overshoot 
(ISamadi et al.ll2010l: IShiode et al.l 120131) . Generally those theoretical models assume Reynolds 
stresses in the convection zone that generate waves with a predominant frequency given by the 
convective turnover frequency and with frequency and wavelength spectra dictated by the assumed 
properties of the turbulence. The amplitudes of the waves are determined by the efficiency with 
which energy is transferred from convection to waves, and is highly dependent on assumptions 






























































-3- 


about the convecti ve-radiative i n terfac e. D espite these s hortco mings, the theoretical spectra of 
IGW predicted by ISamadi et alJ (120lOh and IShiode et alJ (12013h are consistent with the observed 
frequency ranges of variable OB stars. However, their amplitudes vary sig nificantly and appe^ to 
be inconsistent with observations. For example, the amplitudes predicted inIShiode et al.l (120131) are 
more than an order of magnitude lower than those predicted in Samadi et al. J201(i) and even those 
are too small to explain the detection of stochastically excited gravito-inertial modes in HD 51452 
by iNeiner et al.l (12012h. These theoretical tnqdels also assume that the star is non-rotating and it 
was shown bv iMathis. Neiner. & Tran Minhl (120141) that the Coriolis force has a severe effect on the 
stochastic excitation of waves and hence, their surface amplitudes. In reality, the waves in stars are 
mixed gravito-inertial waves and this should be accounted for in theoretical predictions of surface 
wave amplitude^. Beyond their effect on photometric and spectros copic observations, these mixed 
type waves likely lead to substantial angul ar momentum transport (iRoeers et al.ll2013h and may be 


respo nsible for the Be star phenomenon (lAndolll986l: [Rogers et al .1120131: iLee. Neiner. & Mathis 


2014h . which could not be fully captured in earlier ID analytic models. Clearly, observational con¬ 
straints on the detection of IGW, as well as more realistic multi-D predictions of their observational 
signature, would be highly beneficial. 

Here, in order to make comparisons between theoretical predictions of IGW in massive stars 
and observations, we use the velocity spectra from the two-d imensional (20) nu merical simula¬ 
tions of convectively driven waves in a 3 M© star, described in [Rogers et al.l (12013h . While having 
many shortcomings (dimensionality, degree of turbulence), these simulations provide a more real¬ 
istic benchmark for comparison because they self-consistently generate the waves by convection, 
include all non-linearity and rotation. In order to carry out such simulations some sacrifices had to 
be made. Perhaps the most relevant for the current work, are that the simulations used an enhanced 
thermal diffusivity for numerical stability and considered only one mass. In order to compensate 
for the enhanced thermal damping of the waves on their journey to the stellar surface, the waves 
were driven harder by using an (equivalently) enhanced heat flux. While this is not a perfect 
representation of the waves, it is likely the best that can be done with current computational re¬ 
sources. Because of this shortcoming, wave amplitudes produced by these simulations should be 
treated with some caution. Furthermore, the simulations only considered one r nass and it is well 
known that the IG W frequencies depend weakly on the mass of the star (e.g., [Aerts et al.[[2010 : 
Shiode et al.[[2013[l . Repeating the simulations for various masses and evolutionary stages is be¬ 
yond the scope of the current work, given that each 2D simulation takes approximately a hundred 
thousand processor hours for one set of parameters (mass, age, rotation, diffusion coefficients, lu¬ 
minosity). Such a simulation study will be the subject of future research. In order to compensate 


'For ease of presentation in the following comparisons we will use IGW generically to mean convectively driven 
mixed gravito-intertial waves, sometimes referred to as GIW. 
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for the enhanced thermal diffusion, the convective heat flux adopted in [Rogers et alJ (120131) was 
typically 10"^ times higher than that of a 3 M© star, so the predicted spectra of IGW might appro¬ 
priately resemble those of ~ 30 Mq stars since the wave spectrum is predominantly determined by 
properties of the convection, but this is yet to be investigated. In order to account for variations in 
mass with respect to the 3 M© and for the too high convective flux adopted for the currently avail¬ 
able numerical simulations, we apply one multiplicative factor (an unknown parameter) to scale 
the amplitudes of the waves and another one to rescale frequencies for comparisons with real stars, 
as discussed in more detail below. Despite these caveats, we expect that the overall characteristics 
and shape of the predicted frequency spectra should offer a realistic representation of the effect of 
IGW on observables. 


In Fig.lil we show velocity spectra for three cases from [Rogers et al.l ([2013h : a non-rotating 
star, a rigidly rotating star with an initial rotational frequency of 1.100d“\ and a differentially 
rotating star whose core rotates initially with a factor 1.5 times faster than its envelope, the surface 
rotational frequency being 0 .275 d~^ The s e case s correspond with the simulations labelled as Ul, 
U8, and Dll in Table 1 of [Rogers et al.[ ([2013[) . to which we refer for details. Using a radius 


from the prototypical SPB KIC 10526294 with M = 3.25 Mq and R = 2.215 Rq ([Moravveji et al. 


2015h . the two rotating models give initial equatorial rotation velocities of 31 (non-rigid. Dll) and 
123kms“^ (rigid, U8). The velocity spectra were obtained after more than 100 rotational cycles 
(for U8 and Dll), in a layer near the stellar surface at 99% in radius, and the tangential velocity 
is an average over g) e [0,2;r]. In Fig.[Tl it can be seen that the vertical velocity is many orders 
of magnitude smaller than the tangential velocity, as is expected for gravity waves. Fig. [T] also 
clearly illustrates that non-rigid rotation (Dll) inside the star shifts the power to lower frequencies 
compared to the rigidly or non-rotating cases. Given that the few estimates of core-to-envelope 
rotati on of massive stars point to nonrigid rota tion with a factor 2 to 5 faster core than envelope 
(e.g. JAerts et al.[[2003[:[Pamvatnykh et al.[[2004l) . we use Dll for comparisons with observations. 


Beyond the uncertainties in wave amplitudes and frequencies discussed above, it is not en¬ 
tirely clear how to convert velocity spectra, of the type seen in Fig.iU to brightness variations, 
as observed. To address this, we can use the known amplitude ratio of photometric brightness 
variability due to heat-driven g-mode oscillations in SPBs to velocity variations measured using 
spectral line diagnostics in the same stars. The case of convectively driven waves is, of course, 
different as these waves have shorter lifetimes and some are breaking leading to turbulence. On 
the other hand, the heat-driven prograde or retrograde g modes of SPBs correspond with trapped 
internal resonant gravity waves. The 2D simulations produce convectively driven resonant gravity 
waves with amplitudes in the same range as the non-resonant waves. Therefore, in absence of 
any better conversion, we assume that the brightness fluctuation signatures due to velocity varia¬ 
tions from IGW behave similarly to those from heat-driven g modes in terms of amplitude ratios 
of photometric and spectroscopic observables measured at one instance in time. Under this as- 
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sumption, we use the earefully measured amplitude ratios of (A5)/(A<v>) of g modes to eonvert 
tangential velocity spectra to brightnes s spectra. The measured amplitude ratios for the g-modes in 
the single SPB pulsators considered in lOe Cat & AertsI (120021) have values between roughly 1 and 
10mmag/(kms“*), as deduced from the light curves in the B band and from the centroid velocity 
< V > derived from high-resolution high signal-to-noise time-series spectroscopy (see Chapter 6 


m lAerts et alJl201(l for a definition). Below, we use these observational results to scale wave 
amplitudes in the simulations for comparison with observations. 

In the following sections, we first provide observational evidence of the occurence of IGW in 
yumag space photometry of three young unevolved 0-type stars measured with the CoRoT satellite. 
We then predict line-profile variations due to simulated IGW in order to evaluate if IGW can pro¬ 
vide an alternative explan ation of macroturbulent spectral line broa dening, in addition to coherent 
g-mode oscillations (e.g. JAerts et al.ll2009l: ISimon-Diaz et al.ll2010[) . 


3. Photometric signatures: power excess in space photometry of hot massive stars 


The MOST, CoRoT, and Kepler missions revealed a much larger diversity in variability of OB- 
type stars at umag level than a nticipated prior to the era of high-precision space photometry (e.g.. 


Neiner et al. 


201 luAertsIbOlSl for an update). While most of the variability is well understood as 


due to heat-driven stellar oscillations, rotational modulation, magnetic activity, mass loss, various 
sorts of binarity, or a combination of all those processes, the “red noise” powe r excess found in the 
amplitude spectra of three O stars, remains unexplained (IBlomme et al.ll201 lb . 


As is common in the study of heat-driven oscillations, we represent the effect of IGW by 
means of amplitude spectra, where we focus on the ^-averaged tangential velocity (see Figure 1). 
Figure[2] shows the observed amplitude spectra resulting from brightness measurements assem¬ 
bled by the CoRoT mission for three O stars, with surface rotation frequencies of 0.144d“‘ for 
HD 46150 (top panel), 0.249 d“^ for HD 46223 (middle panel), and 0.084 d“^ for HD 46966 (bot 


tom p anel). These spectra are presented in mmag and have been reproduced from iBlomme et al. 


(120111) . who found them to reveal frequencies connected with phenomena of short lifetimes (order 
of hours to days). We compare those observed spectra (grey line in Fig. [21) with the simulated IGW 
spectra shown in Fig.[T]for the non-rigidly rotating case Dll (red line in Fig. [I]). 

We have applied a correction factor of ~0.75 to the frequencies to account for the factor 
~10 higher mass between the observed and simulated stars, following IShiode et al.l (1201 SL their 
Table 1). We realise that the propagation cavity of the IGW changes somewhat as a function of 
mass and age, which may imply that such a simple scaling is not optimal. On the other hand, the 
spectrum of IGW is dominantly determined by the convective flux and less so by the shape of the 
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Briint-Vaisala frequency. For this reason, we consider a simple scaling to be the best approach 
while awaiting future 2D/3D simulations for various stellar masses and evolutionary stages. 

We further applied a factor of 2, 3, and 1.5, respectively to the IGW amplitude spectra (shown 
in red in Fig.O to convert between simulated tangential velocity and brightness, consistent with 
the observed range of amplitude ratios from (AB)/(A < v >) for the heat-driven modes of SPBs. 
The lowest frequencies, below 1 d"\ are not reproduced by the simulations. This could be due to 
too low a stellar mass and/or too low level of differential rotation and/or the too simplified scalings 
we have adopted. However, the overall resemblance is good in terms of frequency structure. This 
interpretation of the “red noise” power excess as due to IGW is currently the only theoretical 
explanation that provides amplitude spectra whose shape is in agreement with those observed in 
the three O stars. 


4. Spectroscopic signatures: macroturbulence in spectral line profiles of hot massive stars 


High resolution time series spectroscopy of confirmed g-mode pulsators show line-profile 
variations (LPVs) of a complex nature b ut such data is only available for a handful of SPBs (e.g., 
De Cat & Aertsll2002l: IPapics et al.ll2012h . In a bsence of systema tic time-resolved LPV studies of 
gravity-mode pulsators among OB-type stars, lAerts et al.l (120091) performed simulations to show 
that the collective effect of numerous heat-driven g-modes leads to LPVs that mimic the effect of 
a macroturbulent velocity (Vmacro)- More recently, it was shown from observed LPVs that even one 
mode o r spot can giye ris e to spectral line broadening in B stars consistent with considerable vmacro 
values (lAerts et al.ll2014t) . These findings indicate that pulsational velocities due to heat-driven 
g-modes could explain the occurrence of macroturbulence along most of the upper main sequence, 
except for the high est-mass core-hydrogen burning O stars which do not show heat-driven modes 
(ISim6n-Dfazll2015b . However, they do show a power excess of frequencies with short lifetime 
(Fig.O, which we demonstrated may be the signature of IGW. It is then reasonable to investigate 
whether IGW can lead to LPVs in these stars. 


Using stellar parameters and frequencies of HD 46150 (IBlomme et al.l 1201 ll: Martins et al. 


2015h and velocity fluctuations given in the Dll simulation (scaled by a factor of 2 as in the red 
curve in th e upper pane l of Fig .[2l) we produced LPVs for its O III 5922A line following the method 
outlined in lAerts et al.l (1201(1 Chapter 6). In doing so, we used the radial and tangential velocity 
fluctuations produced in the simulations (Fig.[T]) as radial and tangential amplitude and we adopted 
a longitudinal and latitudinal dependence according to dipolar prograde spherical harmonics seen 
under an inclination angle of 60°. This procedure produced the LPVs shown in the left panel of 
Fig.[3l Since there are no observed LPVs for these 0-type stars, for comparison we show simulated 
LPVs for the SPB star KIC 7760680, representative of its Mg II4481A line and based upon its 
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identified heat-driven modes detected in the Kepler data. In ord er to generate these LPVs, we have 
taken the frequency, amplitude, and phase values from Table 1 in IPapics et al.l (120151) . The resulting 
LPVs are shown in the right panel of Fig.[3l We find that the velocity fields due to IGW that give 
good agreement with the CoRoT photometric data (Fig.|21) also lead to LPVs that are quite similar 
to those produced by heat-driven gravity modes of KIC 7760680 and that the spectral line wings 
of both stars are predicted to be broadened compared to the case of rotational line broadening 
alone, indicative of macro turbulence. Of course, we have presented just one set of simulated line 
profiles based on particular assumptions about the nature of IGW in terms of amplitude and mode 
type, but these simulati ons are in agreement with the need for macroturbulence to fit the profiles of 
HD 46150, as found by Martins et al.l (120151 Vmacro - 38kms“*). Our simulations also show that 
the anticipated LPVs resulting from the chosen angular velocity dependence and inclination angle 
are tiny in velocity space (v-axis in Fig. [3]) while better visible in relative flux changes (y-axis), i.e. 
the profiles remain fairly symmetric and fluctuate in relative flux level. So actual future detection 
of time-variability in the OIII 5922A line of HD 46150 will require very high-resolution, high 
signal-to-noise time series spectroscopy in order to observe the fluctuations in the broadening of 
the spectral line over time, with respect to a time-averaged profile. 


5. Discussion 


We have shown that the 2D simulations of IGW by [Rogers et al.l (120131) give rise to averaged 
tangential velocity spectra in agreement with measurements of ampitude spectra of flux variations 
in three mid-0 type stars that do not exhibit heat-driven oscillations. Our simulations of LPVs 
based on the velocity spectra due to IGW reveal that such waves can give rise to detectable time- 
dependent line broadening similar to that observed in SPBs. IGW thus offer a viable explanation 
for the occurrence of macroturbulence in the hottest main-sequence stars of spectral type O, which 
are not subject to heat-driven oscillation modes and/or strictly periodic rotational modulation due 
to surface spots. Given the shortcomings of the numerical simulations, it is remarkable that they 
match the observations within factors of two in both frequency and amplitudes. This provides hope 
that numerical simulations may be able to capture the gross properties of waves in stars and that 
more sophisticated simulations will be able to make comparisons more robust. 

Given our results, it is worthwhile to search for signatures of IGW in the residual light curves 
prewhitened for the detected heat-driven mode frequencies for all OB stars observed by the Ke¬ 
pler and CoRoT missions, keeping in mind the recent finding that many of the frequencies found 
in the amplitude spectra of gravity-mode pulsat ors are combination frequencies (po ssibly due to 
nonlinear mode coupling) of heat-driven modes ( Papics et al. 2015 : Kurtz et ^12015 ). However, it 
should be possible to distinguish between those combination frequencies generated by heat-driven 
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g-modes and the signatures of IGW sinee IGW have short lifetimes leading to frequency features 
with broad wings in the Fourier transform, while combinations frequencies due to heat-driven 
modes with infinite lifetimes reveal delta-peaks in the frequency spectra. While we provide this 
interpretation of the nature of macroturbulence in Q stars, an independent s pectro scopic study on 
the nature of macroturbulence has been initiated by ISimon-Dfaz & Herrerd (12014) . 


In the near future , we also plan to ga ther new photometric data with the refurbished Kepler 
mission, baptised K2 (IHowell et al.l 120141) . accompanied by high-resolution high-precision spec¬ 
troscopic time series. Such a combination of data would allow us to test our findings and to provide 
a large enough sample of OB-type stars in all evolutionary phases to test the occurrence of IGW 
across the evolution of the most massive stars. 
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Fig. 1.— Velocity spectra due to IGW for three cases: a non-rotating star (black), a rigidly rotating 
star with initial rotational frequency of 1.1 d“^ (grey), and a differentially rotating star with the core 
rotating 1.5 times faster than the envelope for an initial surface rotational frequency of 0.275 d“^ 
(red), shifted along the y-axis for visibility purposes. Left: tangential velocity; right: vertical 
velocity. 
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Fig. 2.— Measured amplitude speetra (grey) overplotted with the predietions of IGW for ease 
Dll (red) for HD46150 (top), HD46223 (middle), and HD46966 (bottom). The ratios between 
brightness and tangential veloeity variations used are 2, 3, and 1.5, from top to bottom. 
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Fig. 3.— Simulated line profile variations for the OIII 5922A line of HD 46150 assuming IGW 
(left panel) and for the Mg II4481A line of the SPB KIC 7760680, based on its dipole heat-driven 
modes detected and identified in the Kepler data (right panel). The red dashed lines are for rota¬ 
tional broadening alone. 




